Analytical calculation of the transition to complete phase synchronization in coupled 

oscillators 
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Here we present a system of coupled phase oscillators with nearest neighbors coupling, which we 
study for different boundary conditions. We concentrate at the transition to total synchronization. 
We are able to develop exact solutions for the value of the coupling parameter when the system 
becomes completely synchronized, for the case of periodic boundary conditions as well as for an 
open chain with fixed ends. We compare the results with those calculated numerically. 
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I. INTRODUCTION 



The fact that systems of coupled oscillators can de- 
scribe problems in physics, chemistry, biology, neuro- 
science and many other disciplines, is already widely 
accepted by the scientific community. They have been 
used to model diverse phenomena such as: Josephson 
junction arrays, multimode lasers, vortex dynamics in 
fluids, biological information processes, neurodynamics 
[1, 2, 3, 4, 5, J, 0, The coupled oscillators have been 
observed to synchronize themselves (say to a common fre- 
quency value) in a variety of ways, such as total, partial, 
generalized, lag, etc., when the coupling strength varies 
[a, E, ini, 111 111 E 111, [l3- In spite of the diversity 
of dynamics, all these systems synchronize themselves to 
a common frequency, in a tree-like clustering behavior, 
when the coupling strength between these oscillators in- 
creases [11, [13, III IM IM 111 il. Among the coupled 
oscillators systems, is the system of coupled limit cycle s, 
which has attracted interests. The synchronization phe- 
nomena of such system, in spite of the diversity of dy- 
namics can be described using simple models of coupled 
phase equations. 

Although there has been an extensive exploration of 
the dynamical behavior of the coupled limit cycles that 
show the synchronization phenomena, many interesting 
features remain unknown and little has been done ana- 
lytically. Special attention has been given to the com- 
plex synchronization tree of a system of non-chaotic os- 
cillators with nearest neighbors sinusoidal interactions, 
which becomes chaotic when the interaction is turned on 
[24, 25, 26, 27, 28]. This system, a diffusive version of the 
Kuramoto model [1], in spite of its simplicity, possesses 
all the features of phase synchronization of a system of 



chaotic oscillators. 

k = LU., + y [sin(0,+i - e^) + sin(0,_i - Oi)]. (1) 

Here tOi are the natural frequencies, selected randomly 
from a normal Gaussian distribution, K is the coupling 
strength, 9i is the instantaneous phase, 9i is the instan- 
taneous frequency and i — 1,2, . . . , N . Such oscillators 
with nearest neighbors interaction have been seen in the 
literature to describe Josephson junctions, laser arrays 
and phase-locked loops [H, [2^ . These nonidentical os- 
cillators cluster in time averaged frequency, until they 
completely synchronize to a common value of the aver- 
age frequency MMMMM, 0r) ^ (^j), * ^ J, 
where 



= lim — 



^{t)dt. 



Using periodic boundary conditions 
ing the frequencies such that: 



-N — 



N 



(2) 

9i, and seal- 
(3) 



the above system (f ) of N oscillators has a critical cou- 
pling strength K — Kc, where aX K > Kc a complete 
frequency synchronization can be observed and each 9i 
is locked to a fixed value. There is no phase locking for 
K < Kc although the system has clusters of oscillators 
of the same average frequencies. Reducing further the 
coupling strength K , the number of units which cluster 
to the same average frequency decreases, increasing the 
number of clusters, i.e. the number of branches in the 
synchronization tree, until finally all oscillators are inde- 
pendent and acquire their natural frequencies. 

The scaling of frequencies ([3]) in the case of system 
^ limits the synchronization to zero frequency. When 
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this scaling is not used [1^, [13, , the general features 
of the system and the value of will not change, the 
synchronization occurs via the same transition tree but 
shifted to a common frequency value 

1 ^ 

i=l 

as long as we maintain the periodic boundary conditions. 
When we lift them, the system becomes dependent on the 
value of Wo [1^, therefore, we shall keep wq different from 
zero throughout this work. 

In Fig. [T] we show the synchronization tree for a peri- 
odic system with TV = 10 oscillators. Also in Fig. [2] we 
plot the first three largest Lyapunov exponents as a func- 
tion of the coupling strength, K. At and above the criti- 
cal coupling strength, Kc, the largest Lyapunov exponent 
becomes zero with the rest of them being negative, an 
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FIG. 1: (Color online) Synchronization tree for a system of 
= 10 oscillators, with periodic boundary conditions [Kc ~ 
4.70). 
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FIG. 2: (Color online) The first three largest Lyapunov ex- 
ponents for a system with A'^ = 10 oscillators 

indication of synchronization. The initial frequencies for 
Figs. [Hand [Hare chosen as uji = {0.5331, 2.8577, 2.6978, 
2.0773, 2.8257, 1.6062, 0.7788, 1.4680, 1.7622, 2.0190}. 

The analysis of synchronization and the loss of stabil- 
ity becomes particularly difficult since, as we mentioned 



above, it is on long time average that these systems syn- 
chronize. Therefore, for a system (1), one can performs 
analytic calculations only where the variables (phase and 
frequency) become time independent. In order to do 
that, it is necessary to make sure that we are in a zone 
of stationary states so that averages are equal to instan- 
taneous values of the phases and frequencies. Hencefor- 
ward, it becomes clear that we need to arrive at Kc from 
above {K > Kc). 

II. METHOD OF LAGRANGE MULTIPLIERS 

A. Periodic Boundary Conditions 

We shall start with the case of periodic boundary con- 
ditions. At complete synchronization all average fre- 
quencies are equal and all phase differences are constant, 
while its sum equals zero. Since luq is different from zero 
it is necessary to use the phase differences defined as 
01 = &i — O'i-i- Then, the equation ([T]) can be rewritten 
as 

K 

^ivi- Ui^i + y (sin0j_i - 2sin4 + sin^^+i) (5) 

In order to have complete frequency synchronization, the 
above systems of equations should have stable steady 
states. It can be seen from Equation ([5]) that the sta- 
bility of the frequency steady states is independent of 
the initial frequencies as well as the coupling constant K 
provided K > 0. Thus, if wc could establish a criteria 
for the existence of a stable steady states for a minimum 
value of iiT > 0, then simply the minimum value of K is 
equal to the critical coupling strength, Kc, since the sys- 
tem remains synchronized once full synchronization has 
been achieved (for K > Kc). To find this minimum we 
use the method of Lagrange multipliers. 

We should note that the frequencies of the system syn- 
chronize on time average, and after complete frequency 
synchronization all (f)i = 0, and the 0i's are all constants. 
Therefore as long as the calculation is done above Kc, 
we can exchange the variables by their time independent 
values. It is straightforward to see from the steady states 
of (O, where (pi — 0, that 



sin 02 — sin (jjl ' 

and 

sin 0* — sin 02 ^ — ^(sin 0* — sin 0*) = 0, i — 3, i, . . . , N, 

(7) 

where 

a2 = 3(t^o - wi), 

and 

= Oi-i -I- 3(wo - Wi-i), I = 3,4, . . . , A^. (8) 
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In addition, from the definition of (hi, it follows that 



N 



We optimize the fmiction © subject to the conditions ([7]) 



and ([9]) . Then the minimum value of if > is the critical 
value provided that the set {</>*, i = l,2,...,iV} 
form a stable steady state of ([5]). 

For the above purpose, let us define a function 
E{(t)*i,cl)*2, ■ ■ * ,(l>*N,>^o, A3, A4, . . . , Aat) as 



E = 



02 



sin 



— sin ( 



N 
i=3 



sin (jji — sin ( 



r 



ai , . 
■ — sm ( 

02^ ' 



1 ^ 
sin(/)i) +AoE' 



(10) 



where Ai, i = 3,4...,iV are parameters. To have an 
optimum K one should solve 



dE 



dE 



0, i = 1,2,3, ...,iV 



_=0, i = 3,4,...,7V 



(11) 



The above conditions yield the following algebraic equa- 
tion in addition to Eqs. ([7|) and ([9]) 



Recently, Daniels et al. [25| presented a calculation of 
equation (5) to estimate the value of the critical coupling 
Kc based on the results by Strogatz and MiroUo Q . Our 
method, which makes use of the fact that above total syn- 
chronization {K > Kc) all variables: phases and frequen- 
cies, are time independent and the problem of handling 
temporal averages to estimate the value of the critical 
coupling Kc, disappears. In addition, we are able to ob- 
tain a condition for complete synchronization to occur 
[equations dH), ^ and (fl^ ]. 



N 

\ ^ COS ( 



= 0, 



(12) 



Fixed Ends 



We can obtain the steady states (/)* , i — 1,2, . . . , N values 
by solving Eqs. ([7]), Q and (HI]). Hence we solve the 
above set of equations using Newton-Raphson method. 
Assuming a random set of initial values for 0*'s, we look 
for converged and stable values for (j>*'s, then the values 
of if > are estimated using Eq. ©. The application of 
these equations bring the results shown in Table [D 



TABLE I: Critical K values (Kc) for the case periodic bound- 
ary conditions 



N 



Numerical" Analytical' 


4 


2.540775 


2.540771 


5 


2.535587 


2.535583 


6 


3.293280 


3.293245 


7 


4.561624 


4.561588 


10 


4.696059 


4.696048 


15 


2.851727 


2.851651 


25 


3.642530 


3.642538 


50 


9.458391 


9.458378 


100 


12.723401 


12.723417 



"from direct numerical simulation of JTJl 
''by solving Eqs. lO, lO and JH} 



From Table [H we notice that the analytically and nu- 
merically calculated values are in a good agreement. 



We move now to the system where the phases at the 
boundaries are fixed. In this case, with the boundary 
conditions On — cn and = cq, we can write the fol- 
lowing equations: 



-1 = 



(13a) 



K 



= 



(13c) 



where i — 1,2, . . . , N . The equation p3)) can be written 
in terms of the phase differences (f>i — 9i — 0i-i as 

K 

— ivi + — (sin 02 — sin ) , 
K 

= LUi - LUi-i + y (sin(/ii_i - 2sin0i + sm(f>i+i) , 



K 

tiN+i = -^N - y (sin0Ar+i - sin (j)N) 



(14) 



where i = 2, 3, . . . , iV. Then, the constraint in this model 
reads as 



JV-l-l 
i=l 



(15) 



The synchronization diagram is shown in Fig. [3] for 
iV = 10 oscillators with initial frequencies uji = {2.5331, 
2.8577, 2.6978, 2.0773, 2.8257, 1.6062, 2.7788, 1.4680, 
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FIG. 3: (Color online) Synchronization tree for a system of 
A'^ = 10 oscillators, with fixed ends (Ks ~ 5.18 and ^ 
41.13). 



X.IQII^ 2.0190} and the boundary conditions cq = 0.15 
and Cat = 2.15. We notice that the osciUators synchro- 
nize with a frequency just below o^o at a couphng constant 
Ks- Moreover Kg depends on the initial frequencies and 
the values of the phases at the fixed ends, i.e. cq and 
Cat. Increasing the coupling constant, the oscillators re- 



main synchronized in averaged frequency although this 
value of the average frequency decreases as K increases. 
Finally K — Kc, the average frequency collapses to 
zero. The value of K^^ depends on loq as well as on the 
boundary conditions in the phases of the oscillators [1^ . 
From the steady states of Eqs. (fT3)) we find that 



and 



sin — sin 02 



where 



sin 02 — sin (p\ 



— (sin0* — sin0J) = 0, z 

^2 



= 3,4. 



i-l 

^3^t^j, i==2,3,4...,7V. 



(16) 



.,7V, 
(17) 



(18) 



The equation to be minimized by using Lagrange multi- 
pliers is 



E = 



sm ( 



a2 

\ — sin ( 



N 
i=3 



sin( 



sm ( 



— (sm ( 

02 



sm ( 



' N 



[cn - Co) 



(19) 



Upon minimizing (18), we obtain: 



N 

COS (/)* cos 02 



C0S( 



0, 



(20) 



of fixed ends, where it is easy to see that the analytic 
and numeric values of the critical coupling are in a good 
agreement. 



The solution of Eqs. ^ and ^ gives provided 
that we are at stable steady states 0*, i = 1, 2, 3 . . . , A^. 



TABLE II: Critical K values {Kc) for the case fixed ends 
N 

Numerical" Analytical' 

3 30.343908 30.343892 

6 76.993996 76.993767 

10 41.127962 41.127924 

12 68.656599 68.656533 

15 53.059077 53.059054 

25 53.968662 53.967955 

50 115.546391 115.546343 

100 80.218507 80.218133 

"from direct numerical simulation of II13I I 
''by solving Eqs. ([ITJ and (|20]l 

In Table [ll] we show the results obtained for the case 



III. CONCLUSIONS 



In this paper we have performed analytic calculations, 
for a system of coupled oscillators through nearest neigh- 
bor coupling with different boundary conditions, at the 
transition to complete synchronization. These oscillators 
are known to synchronize on time average frequencies, 
forming clusters, increasing the number of oscillators in 
each, until they come to complete synchronization. Un- 
til this is achieved the oscillators remain time dependent, 
chaotic or otherwise, as shown by the Lyapunov exponent 
from Fig. [5J We calculated the value of the coupling 
strength, K^^ at complete synchronization. We found 
that the analytically and numerically calculated values 
are in a good agreement. In addition we obtain condi- 
tions for a complete synchronization to occur. 
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